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We have performed a measurement of the Casimir-Polder force using a magnetically trapped 87 Rb 
Bose-Einstein condensate. By detecting perturbations of the frequency of center-of-mass oscillations 
of the condensate perpendicular to the surface, we are able to detect this force at a distance ~5 um, 
significantly farther than has been previously achieved, and at a precision approaching that needed 
to detect the modification due to thermal radiation. Additionally, this technique provides a limit 
for the presence of non-Newtonian gravity forces in the ~1 fim range. 
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I. INTRODUCTION 

Interest in the Casimir-Polder p| force, the attrac- 
tive QED force between an atom and a surface, and the 
closely related Casimir force, the attractive QED force 
between two surfaces, has blossomed in recent years fol- 
lowing the breakthrough experiments of Sukenik et al. 
and Lamoreaux [E] ■ Additionally, the tremendous exper- 
imental progress in both ultracold atomic systems and 
microelectromechanical systems (MEMS's), has pushed 
both fields towards precise work very close to surfaces — 
regimes where Casimir- type effects become important. 

To the present, experiments have identified the 
crossover in behavior between the van der Waals-London 
and Casimir-Polder regimes, which occurs at an atom- 
surface separation of ~0.1 fim for 87 Rb. Inside of this 
crossover, the van der Waals-London regime, the poten- 
tial scales as 1/d 3 , where d is the distance between the 
atom and surface. Outside of this crossover, the Casimir- 
Polder regime, the potential scales as l/tf 4 . It has been 
predicted that the presence of thermal radiation from 
the surface and surroundings will modify the behavior of 
this force. This modification is predicted to occur at even 
greater atom-surface separations, i.e. ~7 /im at 300 K. 
In this large-separation regime, hereafter referred to as 
the thermal regime, the potential scales as T/d 3 , where 
T is the temperature of the thermal blackbody. 

The experiment of Sukenik et al. |2j was the first to 
clearly measure the crossover from the van der Waals- 
London to the Casimir-Polder regime. A number of ex- 
periments followed that used ultracold atoms to detect 
the presence of the Casimir-Polder force 0, IE IE 0> IE @ ; 
however, none have approached the sensitivity at large 
distances required to detect a crossover to the thermal 
regime. 

Perhaps the most obvious technique for measurements 
of surface forces using ultracold atoms is interferometry 
in which atoms take separate spatial paths. This sort 
of atom interferometry has proved difficult, but some 
groups have now started to make headway in this di- 
rection [IE mill- 



In this experiment, the effects of the surface potential 
on the mechanical motions of a Bose-Einstein conden- 
sate are studied. The attractive surface force distorts the 
trapping potential and thus manifests itself in a number 
of ways. First, the center-of-mass position of the atoms 
changes, but only by ~10 nm for our 228 Hz trapping po- 
tential when the condensate is several microns from the 
surface. Position deviations on this order are well be- 
low our experimental sensitivity; thus, detection of the 
Casimir-Polder force in this manner is unfeasible. Sec- 
ond, the collective oscillation frequencies of the conden- 
sate change. Of these, the center-of-mass oscillation, or 
dipole oscillation, is perhaps the most robust because it 
is very long lived and its frequency is independent of in- 
tracondensate interactions. 

In the most simple approximation, the normalized 
dipole oscillation frequency shift in the x direction [see 
Pig- ID^a)], hereafter referred to as 7 X , can be written as 



> - - -Tr^tr, (i) 

oj x 2uJ x m 

where uj x is the unperturbed trap frequency in the x di- 
rection, lu' x is the perturbed trap frequency, m is the mass 
of 87 Rb, d x is the second partial derivative with respect 
to x, and U* is the potential experienced by the atoms 
due to the surface. Thus the dipole oscillation frequency 
is sensitive to the second derivative of the potential, or 
to force gradients. 

A detailed theoretical analysis of this system was per- 
formed by Antezza et al. [T^j . This analysis includes a 
careful calculation of the Casimir-Polder force from a di- 
electric surface, the modification to the Casimir-Polder 
force in the thermal regime, and the expected normal- 
ized dipole frequency shift j xi taking into account the 
finite width of the condensate and the finite oscillation 
amplitude. The result of this analysis is that for experi- 
mentally plausible conditions, the expected values of j x 
are on the order of 10~ 4 , well within experimental preci- 
sion. 
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FIG. 1: (Color online) (a) Diagram, to scale, illustrating the 
aspect ratio of the condensate and typical oscillation posi- 
tion relative to the surface. The coordinate axis orientation 
and the direction of gravity are also indicated, (b) Typi- 
cal data showing the radial dipole oscillation after expansion 
away from the surface. 



the surface, we use an absorption imaging technique de- 
scribed in [lfl [l7| where we illuminate the atoms with a 
beam perpendicular to the long axis of the condensate. 
This beam impinges on the surface with a slight grazing 
incidence angle of ~2.4° such that when the condensate is 
within ^100 /im of the surface, both a direct absorption 
image and a reflected absorption image of the condensate 
appear. Measuring the distance between these images al- 
lows us to determine the distance between the condensate 
and surface. To calibrate the magnetic field necessary to 
position the condensate a given distance from the surface, 
a series of images are taken where we push the atoms to 
a range of distances ^20-60 ^m from the surface. The 
condensate-surface separations in these images are mea- 
sured and then used for calibration of the magnetic field 
used to push the atoms. 

To allow measurement of surface forces at different sur- 
face locations, the magnetic trap can be moved indepen- 
dently of the surface in the y and z directions. Since the 
condensate only interacts with a ~200 x 10 /mi region 
of the surface, we can translate the trap to measure sur- 
face forces at many different locations on our 5x8 mm 
surfaces. Finally, we can adjust the angle of the z trap 
axis to be parallel with respect to the surfaces. Using 
the surface reflection images, we have verified that the 
deviation from parallel is < 0.25°. 

To excite a condensate dipole oscillation in the x direc- 
tion, we apply an oscillating magnetic field of the form 



II. EXPERIMENT 

We briefly review the apparatus for generating conden- 
sates and measuring surface forces, as a more detailed de- 
scription of the apparatus used to produce the condensate 
can be found in [l4| and the technology and techniques 
for atom-surface measurements are described in detail 
in [TEl Il6[ . At the end of evaporation, nearly pure con- 
densates (the fraction of atoms in the condensate > 0.8) 
of 1.4 x 10 5 magnetically trapped 87 Rb atoms are cre- 
ated in the \F = l,mjf = — 1) ground state. In our Ioffe- 
Pritchard-type magnetic trap, with trapping frequencies 
of 6.4 Hz in the axial direction (z) and 228 Hz in the ra- 
dial directions (x and y), this corresponds to condensate 
Thomas-Fermi radii of 85.9 /im and 2.40 /jm in the axial 
and radial directions, respectively. See Fig. HJa) for the 
coordinate definitions and orientations of the surface and 
condensate in the experiment. 

The surfaces for study are located ~1 mm above (+se 
direction) where evaporation occurs. To position the con- 
densate near the surface, a vertical (x direction) magnetic 
field is applied. This uniform magnetic field acts to dis- 
place the magnetic minimum of the trapping field. By 
applying a carefully controlled field ramp, we are able to 
move the atoms arbitrarily close to the surface without 
exciting mechanical oscillations of the condensate, and 
the condensate can be held there for many seconds. 

To measure the distance between the condensate and 



B a (t)Ke-<*-*> a / Ta cos(w a i), (2) 

where r is the time width of the pulse (10 ms in this ex- 
periment) and to is the time of the peak of the pulse. In 
frequency space, this excitation is centered on the radial 
trap frequency lo x and contains no dc or high-frequency 
components; this prevents excitation of unwanted inter- 
nal condensate modes. Similarly, dipole oscillations can 
be excited in the y and z directions. 

Expansion of the oscillating condensate is accom- 
plished by a microwave adiabatic rapid passage to the 
\F = 2,mp = —2) state, which is antitrap ped , followed 
by ^5 ms of rapid antitrapped expansion |14| . The an- 
titrapped expansion acts to push atoms away from the 
magnetic minimum, and because of gravitational sag, 
the condensate begins the expansion below the magnetic 
minimum, so the condensate is pushed away from the 
surface during expansion. Additionally, the antitrapped 
expansion acts to amplify the radial dipole oscillation 
amplitude by approximately 20-fold, permitting straight- 
forward measurement of the oscillation in expansion. For 
example, see Fig.^b). Finally, the condensate is simul- 
taneously imaged through absorption along both the y 
and z directions, allowing us to monitor the position of 
the condensate in all three directions. 

The typical experiment is performed as follows. First, 
a surface calibration set is taken to determine the mag- 
netic field necessary to position the condensate the de- 
sired distance from the surface. Second, a vertical os- 
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dilation data set is taken at the desired trap-center to 
surface distance d, typically 6-12 /im. Interspersed with 
these data are vertical oscillation data taken at do, the 
distance we use to obtain the normalization frequency u) x . 
Data points and normalization points were randomly al- 
ternated during the course of the data set in order to 
prevent trap frequency drift from affecting our measure- 
ment. For this experiment do = 15 fim. A distance 
of 15 /im is far away enough such that surface forces will 
not affect the frequency; the normalized dipole frequency 
shift from the Casimir-Polder force is less than 1CP 6 at 
this distance. By comparing the frequency measured at 
d to that measured at do, we obtain j x at the particu- 
lar condensate-surface separation d. Last, a second sur- 
face calibration set is taken to determine the condensate- 
surface distance drift over the course of the data set, typ- 
ically <C 1 /^trn. 



III. SYSTEM ATICS 

Rejection of the presence of spurious forces on the con- 
densate caused by surface-based electric and magnetic 
fields is critical to the interpretation of our results. Our 
most powerful test for spurious forces is provided by 
the elongated geometry of our condensate. The mean 
condensate-surface separation for our closest measure- 
ments, ~6 /xm, is significantly smaller than the axial ex- 
tent of the condensate, ^170 /im [see Fig.^a)]. Thus a 
spatially inhomogeneous force, due to localized electric or 
magnetic surface contamination, would likely affect only 
part of the condensate. When oscillating in this inho- 
mogeneous potential, the condensate behaves more like a 
string than a stiff bar and thus will oscillate at different 
radial frequencies along its axial extent. Using a tech- 
nique fully described in |l6j . we analyze images of the 
oscillating condensate and obtain j. x (z) axially along the 
center ^120 /im of the condensate, in addition to 7^ m for 
the center of mass of the condensate. 

To investigate the presence of a spatially inhomoge- 
neous force, we examine 'Jx(z) across the condensate. We 
first define the standard deviation of y x (z) along the ax- 
ial extent of the condensate to be <5 7 . If 5 7 is greater than 
a predefined threshold value 0], then we surmise that 
there is a statistically significant spatially inhomogeneous 
force acting on the condensate and thus move to a new 
surface location. On the other hand, if <$ 7 is less than 
the threshold value, then we take <5 7 to be our system- 
atic limit on spatially inhomogeneous forces experienced 
by the condensate. 

Spatially uniform spurious shifts along the extent of 
the condensate are less likely but must also be accounted 
for. A completely uniform surface charge, or magneti- 
zation, will by symmetry not generate a force. The re- 
maining possible cause of spurious uniform forces is then 
stripes of surface contaminations collinear with the axis 
of the condensate [l^. To test for this possibility, we 
perform measurements at multiple surface locations, as 



well as a series of measurements to test for the presence 
of magnetic and electric fields. 

An atom in an electric field will experience an energy 
shift according to Ue — — (c*o/2)i? 2 , where ao is the 
ground-state dc polarizability and E is the electric field 
magnitude. Thus, to first approximation, the normalized 
frequency shift caused by an electric field can be written 
as 

lx (x -dlU E = ^d 2 x [{E* x f + (E* y f + {Elf} . (3) 

Our goal is then to determine the x dependence of the 
surface electric fields El , E* , and E* ; from these we can 
obtain an estimate of j x . In our previous work [16j , we 
applied a uniform dc external electric field E^ xt , and by 
measuring a change in -f x as E% xt was varied ±100 V/cm, 
we were able to obtain d x E* (x). For the current work, we 
are instead using a technique that allows us to measure 
d x E*{x) as well as d x E*{x) and d x E*{x). 

If rather than apply a dc external electric field we apply 
an ac external electric field E^ xt cos{ojt) , where E ext ^> 
E* , and we invoke Vx£~0 [2(| , then the forces on the 
atoms are 

F x {t) ~ a E° xt cos(u)t)d x E*, (4) 
Fyit) ~ aoE e x xt co S {ujt)d x E* y , (5) 
F z {t) ~ aoE^ xt cos{ojt)d x E* z . (6) 

If uj is set to lo x , ujy, or uj z and d x E* y ^ is nonzero, then 
the oscillating electric field will resonantly drive a dipole 
oscillation. This allows the measurement of very small 
electric forces. By keeping the drive time short, 50-100 
ms in our experiment, compared to the oscillation damp- 
ing rate, typically 1-10 s, and the inverse of the drive 
detuning, typically < 0.5 Hz, then the system reduces 
to that of an undamped, resonantly driven oscillator. In 
this case, the oscillation amplitude linearly grows as 



where Fo is the amplitude of the driving force and A is the 
rate of growth of the amplitude. By positioning the con- 
densate at a given distance from the surface and driving 
a large external electric field (~100 V/cm) at the i trap 
frequency, we are able to detect d x E* , where i represents 
the x , y, or z direction. This measurement is performed 
at multiple distances from the surface, and d x E*{x) is fit 
to a form d x E*{x) = —piCi/ 'x Pi+1 . This corresponds to 
an electric field of the form E* {x) = Cj /x Pi . When fitting 
this data for Ci we vary the exponent pi between 0.20, 
corresponding to a very broad collinear surface charge or 
dipole distribution, and 2.0, corresponding to a line of 
dipoles. Finally, using this range of powers, the power- 
law fit leading to the largest systematic error is used. 
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Although less prevalent, evidence of magnetic surface 
contaminants can be seen on our insulating surface (likely 
remnants from surface polishing). Our trapping potential 
itself is magnetic, so we cannot rely on techniques similar 
to those used for the detection of electric fields. Instead, 
we carefully examine the magnetic trapping potential it- 
self. A spurious magnetic field cannot exclusively modify 
uj x ; rather, the spurious field will manifest itself as a mod- 
ification of multiple trapping frequencies, or anomalous 
center-of-mass displacements. 

The trap frequencies in the three directions can be ex- 
pressed as 



W T — U!„ — 
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(8) 
(9) 



where [Lb is the Bohr magneton, g F is the Lande g fac- 
tor, rj is the linear magnetic field gradient in the x and y 
directions |2l| , Bq is the bias field, and j3 is the magnetic 
field curvature in the z direction \B z (z) = B + (3/2 z 2 ]. 
Adding an additional spurious magnetic field B* and ex- 
panding to first order in B* , the normalized frequency 
shift J x ,y,z induced by the spurious magnetic field can be 
written as 
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Invoking V • B* 
expression for the systematic uncertainty in the normal- 
ized frequency shift caused by a spurious magnetic field, 

8~f Xl as 



8"f X ~ A (S-fy) 2 + 



B 2 {3 2 



(3 2 



(13) 



where we have introduced Sj y and 5j z , which are the 
measured systematic limits on deviations from zero of j y 
and 7 Z as the condensate nears the surface, and 5z cm , 
which is the uncertainty in the movement of the center 
of mass of the condensate in the z direction from its equi- 
librium position. 

This technique can be summarized as follows. If the 
only force on the condensate is the Casimir-Polder force, 
then 7 y , j z , and z cm will not change as the condensate 
is moved near the surface. Therefore, if j y , j z , and z cm 
display no statistically significant deviation, then we have 
verified that magnetic forces from the surface are not 
present at measurable levels. The uncertainty in these 



TABLE I: A summary of our errors in 7^. The la error bars in 
Fig. represent a combination of statistical and systematic 
errors. The relative contributions of the various sources of 
statistical and systematic error were evaluated separately for 
each point. Values for the errors for the worst point and for a 
typical point are tabulated under "Maximum" and "Typical" 
below. 



Error 


Maximum Typical 


source 


(10- 5 ) 


(10- 5 ) 


Statistical 


8.3 


4.0 


Spatial inhomogeneity 


4.6 


2.5 


Uniform magnetic 


2.9 


2.2 


Uniform electric 


4.1 


0.41 


Normalization 


0.18 


0.14 


Total 




5.2 



terms, S^/ yi Sj z , and 5z cmi then allow us to obtain a 
systematic limit on the presence of magnetic forces. 

This technique, where we only consider uncertainties 
of fields, is not applicable in the case of the electrical 
systematics because typically a nonzero electric field is 
detected. Therefore, the perturbation of the measured 
electric field, rather than the uncertainty in the presence 
of an electric field, dominates the electrical systematics. 

Finally, we acknowledge that in certain situations, a 
magnetic field could cause a change in j x and no change 
in 7 y -for instance, if B*/(2B ) = —d y Byjr\. Unlikely 
situations such as this cannot be categorically excluded 
with our current analysis technique; however, by per- 
forming measurements at multiple surface locations and 
verifying that these measurements agree, we can reduce 
the possibility that this sort of unusual cancellation could 
disturb our measurement. 



IV. RESULTS 

The surfaces with which these experiments were per- 
formed were ~ 8x5x2 mm 3 pieces of UV-grade fused sil- 
ica and sapphire polished to ~0.5 nm surface roughness. 
Conducting surfaces would in some way be preferable, 
largely because they are less susceptible to electric fields 
caused by surface charge. Unfortunately, alkali-metal 
atoms, when adsorbed on a conducting surface, gener- 
ate a significant electric dipole field 0] . Preventing any 
atoms from striking the surface during a measurement 
is unfeasible, so is seems that dielectric surfaces, despite 
the possible presence of surface charges, are preferable in 
this case. 

The dielectric surfaces studied include (1) UV-grade 
fused silica prepared by a hydrofluoric acid etch followed 
by UV-ozone cleaning, (2) sapphire prepared with UV- 
ozone cleaning, and (3) UV-grade fused silica cleaned 
with acetone, ultrapure methanol, and de-ionized water. 
Surfaces (1) and (2) displayed forces that were 3-10 times 
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FIG. 2: (Color online) Normalized frequency shift data from 
the fused silica surface (dc dielectric constant = 3.83 (25^. 
Each data point represents a single measurement of y x (no 
data averaging was performed). These data were taken at 
two locations spaced 300 /im apart. Vertical error bars in- 
clude the statistical and systematic errors detailed in Table I. 
Horizontal error bars include statistical uncertainty, surface 
drift, and the effects of the uncertainty in the image mag- 
nification. For this data set the mean oscillation amplitude, 
including the small effects of damping, is 2.06 fim. The mean 
Thomas-Fermi radius in the x direction is 2.40 /im for this 
data. Theory lines, calculated using the theory from Antezza 
et al. dEj, indicate T = K (dashed, black line), T = 300 
K (solid, blue line), and T = 600 K (dotted, red line). Ad- 
ditionally we include the extrapolation of the van der Waals- 
London 1/d 3 potential to these distances (dash-dotted, green 
line) . 



larger than the Casimir-Polder force and displayed signif- 
icant spatial inhomogeneity. Previous studies of a BK7 
surface had led us to believe that magnetic impurities em- 
bedded during the polishing process were a problem, thus 
leading us to try an aggressive surface preparation such 
as a hydrofluoric acid etch. However, with our current 
surfaces, spurious forces appear to be primarily caused 
by spatially inhomogeneous electric surface potentials. 

The fused silica surface (3) displayed the smallest level 
of spurious forces. Nevertheless, even with this surface 
we were forced to study multiple spatial locations in order 
to locate suitable positions for measurements. Suitable 
locations were primarily identified by the criteria that S 7 
be less than a certain predefined threshold value. Ap- 
proximately 40% of the surface regions studied displayed 
spatially inhomogeneous forces. The spatially inhomo- 
geneous forces, identified with the technique previously 
described, displayed peak values ~100% larger than the 
Casimir-Polder force. Spatial variations occurred on ^50 
fim distance scales and displayed ~100% percent varia- 
tions in strength. It is possible that finer structure could 
be present, yet not detectable. 

Once a suitable region was identified, we performed the 
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FIG. 3: (Color online) Current short-range Yukawa-type 
force limits in the 0.1-10 fim range. The limits obtained from 
this experiment are shown by the dashed line. The limits 
labeled Mostepanenko, Kapitulnik, and Price are from |2?1 |. 
|27j . and |28j . respecti vely . The limits labeled Lamoreaux (a) 
and (b) , from |2Sft and |3Q| , respectively, are from two different 
analyses of the Lamoreaux experiment 



experimental procedure previously outlined to measure 
7a; . In addition, a significant amount of data was concur- 
rently taken to put limits on spatially uniform electric 
and magnetic forces. Table I summarizes the limits from 
these systematic measurements. The results of our mea- 
surement of the Casimir-Polder force from surface (3) are 
shown in Fig. 

The first thing to note is that our measurement dis- 
tances are deep within the retarded, or Casimir-Polder, 
regime. This is highlighted by the profound disagreement 
between the measured force and the extrapolation of the 
van der Waals-London force to this distance regime, as 
shown in Fig. We do, on the other hand, see good 
agreement with the predicted Casimir-Polder force from 
our fused silica surface. Unfortunately, we currently do 
not have the experimental resolution to discern between 
the T = K Casimir-Polder force and the T — 300 K 
case. Repeating the measurement at an elevated tem- 
perature, however, appears promising; see Fig. [21 for the 
prediction for T = 600 K. At this temperature, the pre- 
dicted j x is larger than nearly all of our data; thus a 
measurement repeated at this temperature should yield 
a significantly larger signal. Additionally, observation of 
a direct correlation between a change in temperature and 
a corresponding increase the Casimir-Polder force would 
clearly demonstrate the thermal regime of the Casimir- 
Polder force. 

If we assume the T = 300 K Casimir-Polder theory 
is correct, then the data in Fig. can additionally be 
used to put limits on short-range scalar-scalar Yukawa 
potentials of the type 
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U Yu k = - J dV^(l + ae -/ A ), (14) 

where G is the Newtonian constant of gravitation, m is 
the mass of rubidium, p is the density of the attracting 
body (in our case the fused silica substrate), r is the 
distance from the rubidium atom to the volume element 
in the substrate, a and A parametrize the short-range 
Yukawa force, and the volume integral is performed over 
the the fused silica substrate. 

For each value of A, we increase a hypothetical a un- 
til the value of the j x predicted due to the hypothetical 
Yukawa force plus the predicted T — 300 K Casimir- 
Polder force is excluded at the 95% level by the data 
as shifted by worst-case assumptions on the systemat- 
ics (for this analysis uncertainties in spurious magnetic 
forces dominate). The limit obtained in this manner is 
plotted in Fig. [3] with the current experimental limits 
in this region. Redesigning the experiment to optimize 
sensitivity to this signal could permit over an order-of- 
magnitude improvement in the short-range force sensi- 
tivity. This improvement could be accomplished by, for 
instance, using a material with significantly higher den- 
sity or by working over a surface where the condensate 
extends over two materials of different density. 

V. CONCLUSION 

In summary, we have performed a precise measurement 
of the Casimir-Polder force at a significantly larger atom- 



surface separation than has previously been achieved. At 
this large atom-surface separation, effects due to thermal 
blackbody photons become important, and extension of 
these measurements to temperatures ^300 K above room 
temperature should allow a clear detection of this effect. 
Additionally, future experiments performed in nonequi- 
librium thermal conditions, such as holding the surface 
at 600 K with the surroundings at 300 K, are predicted 
to observe significant deviations from the equilibrium 
thermal Casimir-Polder case j3l|. Study of the Casimir- 
Polder force in such nonequilibrium situations will hope- 
fully permit a better understanding of this often nonin- 
tuitive force. 

Finally, this experiment has demonstrated promising 
short-range force sensitivity that could provide limits on 
Yukawa-type forces utilizing a significantly different mea- 
surement type-i.e., atom-bulk vs bulk-bulk. Future work 
with ultracold atoms near surfaces utilizing this tech- 
nique and promising atom-interferometry techniques cur- 
rently being developed, should enable better limits to be 
set in the 1-10 /im regime. 
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